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Abstract 

The exact Green's functions of a Luttinger Liquid (LL) connected to impedances are computed 
at zero and finite temperature. It is also shown that if the resistances are equal to the characteristic 
impedance of the Luttinger liquid then the finite Luttinger liquid connected to resistors is equivalent 
to an infinite Luttinger liquid. Impedance boundary conditions (IBC) include also as a special limit 
the case of open boundary conditions, which are explicitly recovered. Finally bosonization for a LL 
with IBC is proven to hold. 
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I. INTRODUCTION 



Bosonization is one of the standard methods for one dimensional quantum field theories^. 
Discovered independently in condensed matter and high-energy physics it is the main tool 
which has allowed to construct the concept of 'Luttinger liquid' inspired by the physics of 
the Tomonaga and Luttinger Hamiltonians^. The LL is a universality class of ID critical 
systems comprising models as important as the Heisenberg spin chain, the Hubbard model 
or the Calogero-Sutherland model. 

Various boundary conditions (abbreviated as BC throughout the paper) have been con- 
sidered in the past for a LL. The earliest studies focused on the infinite system and the 
finite-size LL with periodic boundary conditions (PBC). The properties in both cases can 
be related through a conformal transformation^. In particular the finite-size properties have 
been very useful in conjunction with numerics for extraction of the LL parameters. Later 
more general boundary conditions were also considered: for instance twisted boundary con- 
ditions (TBC) or open boundary conditions (OBC). TBC led to the discovery of even-odd 
effects in a LL and periodicity of permanent currents with the fiux-. The OBC which allow 
description of broken chains were found to have a dramatic effect on critical exponents lead- 
ing to boundary exponents in addition to bulk ones and bridging the physics of the LL to 
that of boundary conformal field theory^"^. More recently dissipative boundary conditions 
were also introduced for the LL and were used to compute transport properties of the LL: 
they describe the coupling of a LL wire to electrodes; they comprise the so-called 'radiative 
boundary conditions' (which relate time and space derivative of the boson fields^) and a 
'chemical potential matching boundary conditions'- . Other dissipative boundary condi- 
tions include the 'Impedance Boundary Conditions' (IBC) introduced by the author-: they 
consist in a LL connected at its boundaries to two impedances. The IBC can actually be 
shown to encompass both 'radiative boundary conditions' and 'chemical potential matching 
boundary conditions' which constitute special cases of the IBC with boundary impedances 
set at half-a-quantum of resistance h/2e^r- 

We will deal in the present paper with the following issues for a LL with IBC: (1) com- 
puting its exact Green's functions and two-point correlators; this will pave the groundwork 
allowing for (2) extending the usual bosonization technique to such a dissipative system. 

Indeed there is no reason to believe that bosonization of such a system is a valid procedure: 
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bosonization is well established in the non-dissipative situations (infinite LL, finite-size LL 
with PBC, TBC or OBC) but a LL connected to resistors is a dissipative system: plasmons 
have now a finite-lifetime. Actually to the author's knowledge bosonization has not been 
shown to be valid for any LL with dissipative boundary condition. One strategy to show the 
validity of bosonization is to start from fermions and then (by considering the anomalous 
current algebra of density operators) to transform the fermionic Hamiltonian into a bosonic 
onei. Such a course is in our case plagued with difficulties related to dissipation: the density 
eigenmodes of the cavity do not form a neat orthogonal basis of states and do not quantize 
as free bosons. Such problems are actually symptoms of non-hermitian Hamiltonian physics 
and the existence of non-trivial self-energies: this is a recurrent issue for open systems 
which is well-known and has led to recourse to biorthogonal bases of states in such diverse 
contexts as mesoscopic transport^, laser physics (leaking cavities in QED where quantization 
of the gauge field in terms of photons breaks down)^, acoustics^, black hole physics^, etc. 
Biorthogonal bases of states lead however for the bosonization program to unnecessary 
complications. 

Nevertheless it will be shown that bosonization does hold for the model at hand (LL with 
IBC). Instead of following the afore mentioned strategy for bosonization we will find more 
convenient to start from a bosonic theory and then fermionize it. As an interesting side 
result of our proof we will compute the exact Green's functions and correlators of the boson 
Hamiltonian. Another interesting side result with potential applications is the finding that 
with suitably chosen resistances the finite LL with impedance boundary conditions (IBC) is 
equivalent to an infinite LL (namely by identity of their Green's functions): in a mesoscopic 
setting the LL in principle can not be abstracted from its surroundings so that the intrinsic 
properties of an infinite LL are not directly accessible. We show a conceptually simple way 
out which exploits the fact that a LL can be viewed as a quantum transmission line. 

OBC will also be shown to be a special case of IBC corresponding to infinite resistances. 

The paper structure will be as follows: 

- (1) in section [III we introduce a model of a (fermionic) LL connected to resistors through 
boundary conditions (impedance boundary conditions IBC); these boundary conditions are 
then recast equivalently in terms of boundary conditions for chiral bosons. 

- (2) In the next two sections [1111 and [TV] we next focus on the bosonic model and compute 
exactly its Green's functions and correlators (at zero and finite temperature). As side results 
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we obtain the exact finite-frequency conductivity of the LL with IBC and we prove the 
equivalence of an infinite LL to a finite-size LL with IBC with suitably chosen resistances. 

- (3) We prove bosonization for our model in section |Vl 

- (4) We compute the fermion correlation functions in IVIi 

Finally an Appendix (Appendix B) is devoted to the single issue of recovering explicitly 
the fermionic Open BC and the Green's function, with results identical to those published 
in the litterature. 

II. MODEL 

A. Notations and definitions. 

Phase fields: We consider throughout the paper the standard LL Hamiltonian which 
written in terms of the usual phase fields reads: 



In terms of fermionic operators the fermion density operator is related to through: 



(Of course such an identification is not fully warranted at this stage but we will show 
later it does hold even for a LL connected to resistors; for the time being we may view the 
relations as defining abstractly the operator p rather than equating it with the operator 
Similar remarks apply for all the operators defined for fermions such as current, etc.) 

The phase field 9 is defined per: 




where the fields (j) and 11 are canonical conjugates 



[0(x,t); U{y,t)] = i5ix-y). 



n = -d^Q 



and 



[0(a;,t); e{y,t)] = z e{x - y) 



where 9 is the Heaviside step function. 
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Chiral fields: The equations of motion of the phase fields: 



dt(j) — —uK dxQ 

imply: 

(0 ± Ke) = ^- dt {(f) ± KQ) . 
The fields 4>±KQ are therefore chiral and we define chiral phase fields and chiral densities: 



They obey: 0±(x, t) — (t>±{x^ ut). Evidently p — po — p+ + P- 
Current: The particle current density is: 



i = — ^5*0 

'71 



= u{p+- p_) 

where the first line follows from current conservation and the others from the equations of 
motion of the phase fields. 

Chiral chemical potentials: 

We define the following operators /x± (they will prove convenient to define our model): 

. , SH 
P±{x,t) = - — — - 
5p±{x,t) 

where functional differentiation with respect to the particle density p±{x,t) has been 
performed. Physically they correspond to chemical potential operators: their average 
value yields the energy needed to add one particle at position x to the chiral density: 
p± — > p± + 6(x). These chiral chemical potentials correspond to the plasma chiral eigen- 
modes of the Luttinger liquid and not to the left or right moving (bare) electrons. An average 
chemical potential can be defined also as: 
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From their definition it follows that: 



hu 

/^± = YP±' (1) 



where we used the relation: 



H = —- dxpl + pl 

J — L l'2 



Therefore the electrical current = ei 



2Zo \ e e ) 

where we have introduced the characteristic impedance of the LL 

As explained ir>^ the LL Hamiltonian is identical with that of a quantum transmission 
line with a characteristic impedance given as above. 



(2) 



B. ModeL 



Our model consists in a LL connected in series with two impedances. We will use a 
description of the LL connected to reservoirs^ which is the exact implementation of the load 
impedances boundary conditions customary for transmission lines or sound waves in tubes. 

We thus assume the following boundary conditions: 



Zs Ze{-L/2, t) = Vs{t) - /^J^, (3) 

Zs and Zo are interface impedances (at respectively the source and the drain) which 
for simplicity will be assumed to be positive real numbers throughout the paper (in other 
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Figure 1: Impedance boundary conditions: the LL wire is connected to two electrodes at voltages 
Vs and Vd through two impedances. 

words they represent resistors; but more general situations could be discussed with complex 
impedances, which is why we stick in this paper to viewing them as impedances). ie{x,t) is 
the current operator, and source and drain are set at a voltage Vs or Vd (see FigH]). The 
Heisenberg picture is assumed so that we work with time-dependent operators. 

As one can see the boundary conditions are tantamount to assuming Ohm's law at the 
boundaries of the system U = RI: the current is proportional to a voltage drop between 
the reservoir and the LL wire and the proportionality constant is just a resistance. In the 
following the source and drain voltages will be set to zero since we want to compute the 
equilibrium Green's function (in the absence of external voltage). 

It is instructive to recast the (equilibrium) boundary conditions in terms of the chiral 
densities. This yields: 



So that: 



Zsue [p+(-L/2)-p_(-L/2)] 
Znue [p+(L/2)-p_(L/2)] 

pA-L/2) 



hu 



2Ke 
hu 



[p+(-L/2)+p_(-L/2)] 



2Ke 



[p+(L/2)+p„(L/2)]. 



Zs + Zq 



P-i-L/2), 



pAL/2) = |^p,(L/2). 
This introduces reflection coefficients for the density: 



rs 



Zs — Zq 
Zs + Zq 



Zd + Zq 



(4) 



(5) 
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These expressions deserve some comment: they are just what one would expect for a classical 
transmission line connected to load and drain impedances. The (classical) equations of 
motion are indeed also valid at the quantum level since the LL Hamiltonian is quadratic so 
we might have anticipated any linear relation to carry on. The basic physics of the boundary 
conditions considered in this paper are therefore those of standing waves in a transmission 
line produced by reflections at the boundaries due to impedance mismatch. 

The reflection coefficients for the phase fields can also be derived; from eq.dH) it follows: 

40+(-L/2,t) = r5 9,.0_(-L/2,t), (6) 

40-(V2,t) = rD9,0+(L/2,t). 

The IBC can then be rewritten as conditions on the non-chiral phase field: 

-ZsUtct>{-L/2,t) = Zod,<P{-L/2,t), (7) 

Zn-dt(t>i-L/2,t) = ZodM-L/2,t). 
u 

As an aside we note that open boundary conditions (OBC) are recovered by setting 
r5 = r£) = 1 (see Appendix B). 



III. GREEN'S FUNCTION OF THE PHASE FIELD. 
A. Green's function of the phase field (p. 

Having derived boundary conditions for the boson fields we now forget the underlying 
fermions and focus on the following model: 

hu [^'^ , 1 _ .._2 



hn C^l^ 1 

2 J-Ln 



L/2 

where as before the fields (p and 11 are canonical conjugates and: 

(f)±KQ 



2 

The problem we will tackle is the following: solving for the retarded Green's function of 
the bosonic Hamiltonian subjected to the boundary conditions for the fields: 

d,M-L/2,t) = rsd,(f)4-L/2,t), 
d,(P_{L/2,t) = rDd,(P+{L/2,t). 
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These boundary conditions introduce dissipation in the problem. 

The Green's function can be conveniently divided into four chiral components: Gr — 
Gp- + 0]^- + G^~ + G^^ where G%'^{x,t;y,t') = -ie{t-t') ([0±(a;, t), 0±(y, t')]) for the re- 
tarded Green's function and appropriate definitions for the advanced Green's function. Since 
we will also need the chiral Green's functions for the bosonization proof instead of directly 
computing the full Green's function we will first compute the chiral Green's functions. 

Using the equal-time commutation relations for the chiral fields, 



[0±(x,i),0±(y,t)] 



{x,t),(l)-{y,t)] = 



iK 

±^ sgn{x -y), 
iK 



one gets the equations of motion for the chiral Green's functions as: 



dx u dt 



K 



G^^{x,t;y,0)^—sgn{x-y) 5{t), 



Au 



and: 



d d' 
dx u dt 



G^^{x,t;y,0)^-^6{t). 



The impedance boundary conditions imply: 



d.G++ = rsd,G-^+ , X = -L/2, 

9xG%~ = rsd^G'^' , x = -L/2, 

dxGjf~ = rnd^G'j^' , x = L/2, 

dxG-j+ = rod^Gp- , x = L/2. 



After Fourier transforming according to: 



m=i ^/Me" 



■ILLlt 



and defining k — uj/u the equations of motion imply the following forms: 
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Likewise: 

K 



Q+- ^ ce^*=^ + 



G-+ = de-'''''-^. 
The boundary conditions imply the following relations: 

= —ik rs e^'^ d 
a> = -ikr^'' e-'^P d 
b< = ikrs^e-''^c 

— ik td e"^ c. 

where we have defined a phase </? corresponding to the phase accumulated along the wire 
the plasma wave: 

(/? = kL. 

Using now the discontinuity of the derivatives of G^^ and G yields: 

2u 
2u 



Finally: 



^ sgn{x - y), 



Aicu 

Gt ix.y.uj) =- —7 ^ ^e'"^^+^^ + -; — , 

Gfj (x.y.u}) =- ^r^e 



2ia; (1 — rgrue*^"^) 4ia; 
Gathering all terms the full Green's function is thus: 



10 



(where </? = kL). As a function of time this yields: 
Giiix, t; y, 0) = -f En=o(^5rz.)" {e {t - + rsmO (t - feU^^) 

-rse(t- (^-+^)^+-+y ^ -rDe(t- (^-+^)^—y ' 



For t < 0, Gr = as it should be. 

Technical note: the studious reader interested in deriving directly the advanced Green's 
function from the boundary conditions should take note that they must be modified. The 
reason is simple: a retarded Green's function which describes outgoing waves (away from 
the system) are clearly compatible with dissipative BC; but an advanced Green's function 
describes incoming waves. So we need BC invariant under time-reversal: to render the time 
derivative compatible with time-reversal we multiply it by sgn{t) 

- sgn{t) Zs dt(t>{-L/2, t) = Z^d,(t>{-L/2, t), (9) 
sgn{t) Zd dt(f){L/2, t) = Zod^(f){L/2, t). 

In terms of the chiral fields this implies: 

d,M-L/2,t) = rT"-'^ dA-{-L/2,t), 
d,<l>4L/2,t) = rT^'^d^MLM- 

One can check that the BC are now compatible with the usual relation 

Ga{x, t] y, 0) = GR{y, -t; x, 0). 

Note also that the two-point correlators mix advanced and retarded Green's functions: 
therefore they will obey these modified boundary conditions as can be readily checked. 



B. Discussion. 



1. Interpretation. 

The interpretation of the Green's function is quite straightforward: to propagate from 
one point to the other there are four kinds of basic trajectories (see Figure 2), (1) one can go 
straight from the starting point to the arrival point, or (2-3) go after bouncing against one of 
the two boundaries, and (4) lastly go after bouncing two times against different boundaries. 
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Figure 2: Interpretation of the Green's function. 

These basic trajectories must then be convoluted by round trips along the whole loop (of 
length 2L) which yield the overall factor (1 — rsr^e^'^'^) ^ (where (p = ^L) in the frequency 
domain expression in eq.(l8|). 

The main difference for the chiral propagators is the appearance of zero modes. The other 
terms have straightforward interpretations: as before they correspond to straight trajectories 
from ?/ to X or to propagation with bouncing at either or both of the boundaries. The 9{x—y) 
or 6{y — x) come from the fact that chiral propagation prevents some trajectories depending 
on the respective positions of x and y. 

We have for now only discussed the retarded Green's function but of course the advanced 
Green's function is simply related to the retarded one through: 
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Ga{x, t; y, 0) = GR{y, -t; x, 0). 

So there is no more work to do. The causal Green's function can also be computed: for 
instance from the two point correlator (which is discussed later). 

We also observe that the retarded and advanced Green's functions do not depend on tem- 
perature since the equation of motion obeyed by the Green's functions has no temperature 
dependence. We have therefore found the phase field Green's functions at all temperatures. 

2. Poles and relation to open boundary LL. 

The open boundary conditions (OBC) are a limiting case of the IBC considered in this 
paper. Hard-walls at the boundaries can be reproduced by considering infinite resistances 
Zg and Zd this then implies that the reflection coefficients are equal to unity rg — td — 
It can then be checked that our expression for the Green's function reduces to that found for 
the LL with open boundary conditions. We refer to Appendix B. Open boundary conditions 
are actually closely related to the ones discussed in this paper: indeed the excitations in 
the LL with IBC are those of the LL with OBC albeit with a finite-lifetime. While in the 
infinite LL or the PBC LL one has travelling waves these excitations are just the standing 
waves expected from a system enclosed within boundaries. 

Indeed the poles are simply: 

uJn = -j^ + M'f^srD)- (10) 

For the OBC there are nodes at the boundaries; the resonances are therefore such that the 
length of the wire is L — nX/2 which leads to cun — Within the simple model of boundary 
resistances (real positive impedances) the lifetime r = 2L/ {u \n{rsrD)) is independent of 
the index mode: the level broadening is constant for each standing wave plasma mode. 
But more complicated situations can be considered: if we assume frequency dependent 
complex impedances the reflection coefficients then acquire a frequency dependence. This 
however does not affect the validity of the expressions in the frequency domain just derived 
for the Green's function. However the structure of the poles will not be quite as simple as 
that described above, since the poles are now determined by: 
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1 - rs(cj)r^(w)e*2W" = 0. 

How can we probe these poles? One of the simplest way is through conductivity or 
conductance measurements. These poles will show up as resonances. Indeed as shown in 
Appendix A the conductivity is simply: 

g 

a{x,y,uj) = -ujGR{x,y,uj). 
vr 

The conductance (which is a matrix in this context of a gated wire connected to two 
electrodes) was computed elsewhere^. 

3. Impedance matching. 

Let us recall the basic physics of transmission lines: an ideal transmission line or LC-line 
(e.g. a coaxial cable) has an energy per unit length £{x) = |£j^(x) + ^p'^{x) where £ 
and C are respectively an inductance and a capacitance per unit length. For an infinite 
transmission line the eigenmodes are traveling waves (plasmons) with velocity u = 1/ \J LC. 
For a finite transmission line connected to two resistors at both boundaries one observes 
reflections: in general any pulse injected in the transmission line is reflected which leads 
to energy losses. In order to minimize losses electrical engineers take advantage of the 
phenomenon of 'impedance matching': if the resistors have identical resistances equal to Zq = 
(the characteristic impedance of the transmission line) then the reflection coefficients 
vanish so that no reflections can occur in the combined system of two resistors+transmission 
line, which becomes effectively lossless. The flnite transmission line has become effectively 
equivalent to an infinite transmission line. This is the origin of the normalized characteristic 
impedance of coaxial cables. 

What is the relation to the Luttinger liquid ? A LL is actually a quantum transmission 
line. Indeed its Hamitonian density is just that of a quantum LC-line since: 

f , hu n huK .9 

which rewritten in terms of the charge density and the charge current pe = ep and je = euKj 
becomes: 
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with: 



C 



h 



C 



2Ke^ 



2uKe^ 



hu 



So it is only natural to inquire whether the 'impedance matching' physics is still valid at 
the quantum level. Quite remarkably it is. We prove the following theorem: 

Theorem: the physics of a finite length Luttinger liquid connected to two resistors 
having resistances equal to the characteristic impedance Zs = = Z^ is equivalent to that 
of an infinite LL (for any observable defined on the length of the LL). 

Proof: the proof follows from identity of the Green's functions. Indeed Zq = Zs = Z^ 
implies: rs = td = 0. The expressions of the one-body Green's functions we have computed 
then trivially reduce to those of an infinite LL. The N-body Green's functions are therefore 
also equal since by Wick's theorem they reduce to a product of single particle Green's 
functions. 

Here in the quantum case the cancellation of the refiection coefficients leads to the 
same physics as in the classical case. Note that the result is still valid if we assume fre- 
quency dependent refiection coefficients: if one checks carefully our derivation of the re- 
tarded Green's function in section [1111 one will notice that the expressions do not require 
frequency-independent coefficients. 

Such a result might be useful in any situation where the intrinsic properties of a LL (infi- 
nite system) are needed: for transport experiments on quantum wires or carbon nanotubes 
the coupling to the leads (source and drain) unavoidably modifies the pure physics of the LL. 
Should one be able to tune the impedance of the leads, one might be able to disenfranchise 
oneself from the interfering effects of the leads. 

In general however the leads will couple to the LL not only with a resistive component but 
also a capacitive (or even an inductive) component. To achieve perfect impedance matching 
would mean to be able for all frequency to adjust Zsiuj) = Zq = Zciuj). For all practical 
purposes depending on the phenomenon one wishes to observe it will be enough to match 
impedances on a finite window. 
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Figure 3: Impedance matching. 

Note also that the possibility to match impedances rests crucially on an independent 
measurement of the characteristic impedance Zq. This can be achieved through a variety 
of means: for instance through the finite-frequency conductance as explained by the author 
in^ or through tunneling experiments. 

However it might not be completely necessary to measure the characteristic impedance: 
one might think of time-resolved experiments where one sends a charged pulse in the wire; if 
one detects reflected pulses this means that impedances are not matched (actually such tests 
are routinely used by telecom operators on their transmission lines to find broken lines). 

IV. TWO-POINT CORRELATORS OF THE PHASE FIELDS. 

We derive in this sections several correlators at zero and finite-temperature. They will 
be needed for the proof of bosonization but are also interesting in themselves since they are 
the building blocks for the vertex operators correlators. All observables of interest can then 
be calculated: tunneling density of states, current noise, etc. 
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A. Correlators of the fields (f) and G. 

It is a simple matter to extract two-point correlation functions; using the definitions of 
the Green's functions one has indeed at zero temperature: 

{(f){x, u;)(f){y, -u;)) = i e{uj) [Gr{x, y, u;) - Ga{x, y, u;)] . 

Tedious but uneventful calculations finally yield the following correlator at zero temper- 
ature : 



F{x,y,t)^ {<f>{x,t)<P{y,0)) 



(11) 



S + i{t + 

u 



+ 



,x-y.2 



u 



+ f.T.:=o{rsrnTrs In J] 

£=±: 

+ E^o {rsmT ro In J] 



5 + <i + e — —) 



u 



e=±l 



5 + i{t + e — '—) 



u 



Extension to finite-temperature is done by observing that the retarded and advanced 
Green's function are temperature independent. Let us consider then the commutator: 



M(x,i;y,o) = mx,t)- 4>{yM) ■ 

- {ct>{x,t)cj>{y,Q)-cj>{y,Q)ct>{x,t)) 

By fluctuation-dissipation (Lehmann's spectral decomposition): 

M{x,y,uj) = (l-e-^'^)((/.(a;,a;)0(|/,-a;))^ 
= {e^- -l){(t>{y,-uj)(t>{x,uj))^ 

at temperature T. At zero temperature this implies: 

M{x, y, uj) = e{uj) (0(a;, uj) <p{y, -uj))^ - e{-uj) {(p{y, -u) uj))^ 

where the subscript denotes zero temperature. Since M is actually temperature independent, 
comparison of the two expressions yields: 

1 



F{x, y, u, T) = {(t){x, u) (t){y, -u))^ 



1 - e-P^ 



[e{u) F(x, y,u;,T^O)- 9{-u) F{y, x,u;,T ^ 
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so that the finite-temperature correction is: 



F{x, y, T)-F{x, y,u;,T^O)^ 9{u;) -^^J—^ F{x, y, u;, T = 0)-e{-u;) ^ F{y, x, to, T 

To regulate the previous correlator one defines: F_+(x,y,t) = {(j){x,t) (j){y,0)) — 
'^{cj>^{x,t) + cl>^{y,Q)). 

The finite-temperature correlator acquires then a correction from its zero temperature 
expression AF = -F_+(x, y, t, T) - F^+(x, y,t,T ^0): 



K 

4:71 



+00 



^ (rsrc)'"' In 



\T [1 + if3-^2nL)](- 



|r [1 + ip-^{t + ^-^+'"^ )] r [1 + i/3-^{t + ^-"+'"^ 



K 



+00 



K 



n=0 





r 


ri + i/3-l(2-+(2n+l)Lj 


2 


r 




2 


r 




2 


r 




2 



+^ (^^s^d)" In 



n=0 





r 




2 


r 




2 


r 




2 


r 




2 



where we have made use of the product expansion of the Gamma function F and P is ^. 
This can also be rexpressed as: 



AF = — 



K 
An 



K 



n=—oo 
+00 



[shcnp-^2nLY 



n=0 

+00 



+ ^ ^ (rsri?)" ro In 



n=0 









shcTT 








+1)L^1 













whare we have defined the function shcX — sinhx/x. 

The correlator for the other phase field © can be found from the previous expressions 
by the operations: {rs,r£,,K) — {—rs,—r£),l/K). The change of sign for the reflection 
coefficients follows from the fact that the chiral components of 6 are 6± — ±(j)±/ K so that 
the refiection coefficients for acquire a relative minus sign with respect to those for (j). 
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B. Chiral correlators. 



The chiral correlators will be useful for computing the cross-correlators of the phase fields. 
We use again: 

{(t)±{x,uj)(t)±{y,-uj)) = ie{uj) [G'^{x,y,uj)-G'^{x,y,uj)\. 

(= -2e{u)^mGi^{x,y,u)). 

Using the expressions of the chiral retarded propagators and since G^{x,y,uj) — 
G'^{y,x, —uj) and G'^{x,y,uj) — G'^{y,x, —u) one finds finally: 



K 



5 + i[t 



2nL ± {x — y) 



u 



K 



+00 



0)) = - J] {TsTnT \ rs In 



n=0 



S + i[t- 



x + y+ (2n+ 1)L 



u 



+rD In 



5 + i it- 



K 



x + y- (2n+ 1)L 



u 



.K .K ts-td 
I h I — 



n=0 



5 + i\t + 



1 - TsTd 

x + y+{2n + l)L 



u 



+rD In 



5 + i{t + 



x + y- (2n+ 1)L 



u 



.K .K rs-TD 
+ 1^ — 1- 



1 - rsTD 



Note the presence of non-trivial zero-mode terms for the cross-correlators of the chiral 
fields. 

For later use in the bosonization proof it will be useful to consider these expressions 
for X = 1. To avoid confusion we define the chiral fields for X = 1 as (pR and ipL\ their 

correlators are therefore as above with K = 1: 

+00 

47r 



^ +00 

(<^R/L(a;,t)93R/L(y,0)) = -— ^ (rsri^)'"' In 



5 + i[t- 



2nL ± {x — y) 



u 



{ipR{x, t)(pL{y, 0)) = ^ X] i^sro f < rs In 

n=0 ^ 



6 + i{t- 



x + y + {2n + 1)L 



+rD In 



^ +00 . 
{iPl{x, t)(pR{y, 0)) = 2- X] ('^sTdT \ rs In 



X + ij - {2n + 1)L 



u 



u 



n=0 



+rD In 



S + i{t + 



x + y- {2n + l)L 



8 8 1-rsrD 
x + y + {2n + l)L 



u 



u 



8 8 1 - rsTD 
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C. Cross correlators of (j) and @. 



Since: 



e 



+ 

K 



one has: 



(06) = K-^ {(0+0+) - (0_0_) - (0+0_) + (0_0+)} 
(0 0) = {(0+0+) - (0_0_) + (0+0_) - (0_0+)} 



Finally: 



(0(a;,i)e(y,O)) 



n=— oo V V u I 

1 +00 r5 + ,^_£±I±(^!!±l)^j 



^ +0O 



a;+j/-(2;) + l)L 



_^ 2, _ Z Ts-Td 



4 4 1 - r^ro 



{e{x,my,Q)) = -T- E (r5r^)'"'ln 



2n L-\-y—:r 



n=— oo ^, V w 

+0O f^ + .ft-^i^i^^^) 
+ ^ y (rsmT rsln { ) f 



47r 

n=0 



47r 

n=0 



1 + 1 ^S-^D 



4 4 1 — r^r^:) 

At finite temperature the regularized correlators Fi — (0(x, t)Q{y, 0)) — | (0(x, 0)©(x, 0)) — 
i (0(y, 0)e(y, 0)) and F, = (e(x, i)0(y, 0)) - § (e(x, O)0(x, 0)) - § (e(y, O)0(y, 0)) get the 
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corrections: 



and: 



ra=— oo \ V u I 



n=—oo \ ^ 1^ \ u I 



n=— oo 

-j^ +00 



47r 

n=0 



1 [sh^-Kii-^ |^^_ a.+,-(2n+l)L \ 



47r 

n=0 



where as before we have defined shjyx) = sinhx/x. 



V. BOSONIZATION. 



It is far from obvious that bosonization works with the boundary conditions considered 
in this paper. 

We now show it does or more precisely we show that: 
The bosonic field theory on a finite length with boundary conditions 

rL/2 ^ ^ 

J-L/2 

d,ipR{-L/2) = rsd,ipL{-L/2) 
d^ipL{L/2) = rDd^(pR{L/2) 

with : 

[pR/L{x,t),^R/L{y,t)] = ±-sgn{x-y) 
[(fiR{x,t),(pL{y,t)] = - - 
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is equivalent to the fermionic field theory: 

rL/2 

H = dx -iR+d^R + iL+d^L 

Pr{-L/2) = rsPL{-L/2) 
Pl{L/2) = rDPR{L/2) 

where R and L are Fermi fields obeying the usual anticommutation rules and: 

PR {x,t) = : R+{x,t)R{x,t) : 
PL {x,t) = : L^{x,t)L{x,t) : 



provided we identify: 

R{x,t) = exp2i\fKipR 
v27r5 

Lix.t) = exp —2i\fKU)L 

V2^ 

where S is a short- distance cut-off. 

Note that we have normalized the Fermi velocity to u — 1. 
We proceed in three steps: 

(1) we check that the vertex operators obey fermionic commutation relations; 

(2) by point-splitting we show that the usual relations still hold: 

: PR {x,t) := -^d^LpR 



: PL {x, t) := -^d^ipL 



(3) and also by point splitting that: 

■■PR {x,t) :2:=^ ■.R+d.R: 
■■PL {x,t) ■.^■.^ ■.L+d.L: 

from which the proof is trivially completed by substitution. 
Step 1: 

For R{x, t) — -^=^exp2i^i^R the relation i); it!(y, i)} = is automatically fulfilled 

given the commutation relation of the chiral field. 
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Indeed using the identity : 
it follows that 

R{x,t)R{y,t) = i?(|/,t)i?(x,t)e-^^<['^«(^'*)'<^«(2/,t)]> 

= -R{y,t)R{x,t) 

Similarly {R{x,t); L{y,t)} = using [(pR{x,t),(pL{y,t)] = -i/4. 

This step is trivial but it remains to show that {R{x , t) , R'^ {y , t)} — 5{x — y) which is 
less straightforward. 

We use the expressions derived for the chiral correlators (inserting K = 1 and u = 1) 
(note that to make clear that the fields are taken at = 1 we use the notations (pR and ipL 
instead of (j)± in the whole section): 

{ipR{x,t)ipR{y,Q))^-— {TsTDt^\n[5 + i{t-2nL±{x-y))\. 

n=— oo 

This entails: 
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{R{x,t),R+{y,t)} = ^ : exp2i^ {^R{x,t) - ipR{y,t)) : 
X {exp An {^r{x, t)(pR{y, t) - ^ {^r{x, t) + t))) + 
+ exp Att {(pR{y, t)ipR{x, t) - | {(fHx, t) + <^\{y, t)))] 
= 2^ : exp2iV7r((^R(x,t) - : 

X {exp - E:r_oo {rsruP In (^±^gSF^) + 
+ exp - E:r_oo {rsroP In (^^^±SS±F^) } 
= 2^ : ea;p2i0r((^ij(x,i) - : 

X {exp -ln(^^) exp - En^o (^^rz.)'"' In (^^) 
+ exp -ln(^±^) exp -E„^o(^5r..)l"l In (2^)} 
= ^ : exp 2i t) - My, t)) : exp - E.^o (rsroT In (^^^) 

X [ ^ + 1 ] 

5-i{x-y) S+i{x-y) J 

= : ea:p2zV^((^«(a:,t) : exp - E„^o (^^r^^)" In - y) 



Step two: 



R'^(x)R(x) — lim - — - exp —2^^/7^(fiR(x + s) exp2i^/7:ipR(x) 

^ ^ini : exp -2^0? (99^(0; + e) - <^ij(a;)) : 
£->o ztto 

exp 47r + e)ipR{x) - ^ + e) + ip%{x)) 

+00 



n=— 00 ^ 

= !^-^o2^ [1 - 2zV^s9.^«(x)] ^ exp-;^(r5r^)l In J 

In the last line the regulator can safely be put to zero in the exponential (except for the 
n = term) since the 2nL term is finite. 5 being the short-distance cut-off (the inverse of 
the bandwidth) is smaller than any distance, so smaller than e. The limit 5 — > is therefore 
taken before the £ — > : 
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R+{x)R{x) = lim — i— [1 - 2iV^ed,Mx)] exp - V (r^r^)'"' In f '^^^f 

and we get after normal ordering (discarding a c-number piece) the expected result. So 

1 

PR {x,t) = : R'^{x,t)R{x,t) := —=d^(pii 
PL {x,t) = : L^{x,t)L{x,t) := —=d^ifL 



Step three: 

We will need the correlators of the vertex operators. 
{R+{x)R{y)) = ^ (: exp -2i^/^ {(fR{x) - (finiy)) :> exp Att LpR{x)^R{y) - ^ + <^^(y)) 



1 / x|n| 1 [ S + i{2nL + x- y) \ 



and: 



(i?(x)i?+(y)> = 2^ (: exp 2iV7r - (^ij(y)) :> exp 47r (^(fR{x)(fR{y) - ^ (<^r(x) + ^^(y)) 

n=— oo ^ ' 

We now compute the square of the normal ordered chiral density. This is of course a 
singular operator which needs itself to be normal ordered. By point splitting and then 
Wick's theorem: 

: R+{x)R{x) :^ = lim : R+{x + e)R{x + e) :: R+{x)R{x) : 

= lim : R+{x + e)R{x + e)R+{x)R{x) : + + e)R{x)) : R{x + e)R+{x) : 

+ R+{x + e)R{x) : {R{x + e)R+{x)) + {R+{x + e)R{x)) {R{x + e)R+{x)) 

= lim—— exp - > [rsTD) In — — 

n=— oo ^ ' 

{:R{x + e)R'^{x) : +R^{x + e)R{x) :} 
+ c — number. 

As before in the exponential only the n = is really singular in the limit 5 — > 0; it cancels 
with the prefactor so that after normal ordering: 
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:: R^ix)Rix) -? := lim exp — (r^r In 

{■.R{x + e)R^{x) : +R^{x + €)R{x) :} 
= lim {■.R(x + e)R^(x) : +R^(x + e)R(x) :] 

= — : R^dxR '■ 

TT 

Finally: using the relations proven in step two the Hamiltonian can be rewritten in terms 
of the currents of the fermion vertex operators: 

rL/2 

H = n Ldx : pI + pI : 

J -L/2 

Then using the relation derived in step three: 

■.pR :2=^ -.R+d^R: 

■.pL {x,t):^=^ -.L+d^L: 

implies immediately: 

rL/2 

H = dx - iR+d^R + zL+S^L : 

J -L/2 

It remains to prove that the boundary conditions for the boson theory translate into the 
quoted boundary conditions for the free fermion theory: but this has already been shown in 
section [Til of the paper (see eq. (l3|4f6l) ). 

Switching on interactions. What about interactions? Having shown the relation 
between the free boson and the Dirac fermions, interactions can now be switched on for the 
fermions. But since we have proven that the customary dictionary of correspondence still 
holds (vertex operator, currents ) it is clear that for the LL with boundary conditions the 
transcription of fermion interactions will go as in the standard LL. For example if we add 
the interaction: 

rL/2 

V = dxg2 pRpL + 94 {pR + PL^) 

J -L/2 

the full Hamiltonian can be rewritten in terms of the phase fields as: 

1 f^l'^ 

H=- dxil + g, + 2g2){dAf + {1 + 9a- 292) 

^ J -L/2 
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f2nL + e\ 
V 2nL J 



which takes the standard form: 

H=n dx-{d^<pf + KI]^ 

^ J-L/2 

with 

K = J' + 3.-2g. 

V 1+^4 + 2^2 

u = sJ{l + g^f-Agl 

(we have normalized the Fermi velocity fi;' to 1 and also h = 1 m. all this section). 

If we look closely at the bosonization proof we see that what does the trick is always 
the n = term in the sum (this term corresponds to reflectionless propagation in the 
propagator): this is sensible since it corresponds to what one would have without reflections 
as in the infinite system. Remarkably this shows that the ultralocal structure of the current 
algebra is unaffected by the dissipative boundaries. 

What about spin? Extension to spinful LL is done in the same manner as with the usual 
LL. We consider a second copy of the LL with IBC and add a spin index a. Given the fact 
that the boson fields (pn^ and ipia for different spins commute, as usual it suffices to add 
Majorana fermions rj^j to enforce equal time commutation relations for the fields for different 
spin&i' so that: 

and 

{i?^(x,t),i?+,(?/,t)} = 5^^> S{x-y) 
if {Va, Va'} = 26aa' aud r]+ =1]^. 



VI. FERMION CORRELATORS. 



The fermion operators are given by the usual relation with the phase fields : 

R{x,t) = -^=expiy/7i{(f) + Q) , 

L{x,t) = expiy/7T{-(p + Q) , 

up to unessential phases corresponding to shifts of the chemical potential (see Appendix B 
and the discussion regarding open boundary conditions). 
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The two-point correlators are given by: 



{R{x, t)R+{y, 0)> = 2^ expTT t)<f>{y, 0)) + {e{x, t)e{y, 0)) + {<P{x, t)Q{y, 0)) + 

+ (e(x, my, 0)) - i t) + 0)) + {e2(x, t) + e^y, O)) + 

+ {<P{x, t)e{x, t) + <P{y, 0)9(1/, 0)) + (e(a;, t)<P{x, t) + e(y, O)0(|/, 0))]} 

and a similar relation for the left fermion. It then suffices to insert the expressions derived 
in the previous section. 

For instance at zero temperature: 



{Rix,t)R+{y,0)) 



hoo 

2TrS J- in=— oo 



{5+i2nLy 



[S+i{ut+2nL)f+{x-y)'2 



K+K~ 



-(rsro)l"l 



X 



n 



oo 
n=— oo 



S+i{ut+2nL+x—y) 
d+i(ut+2nL+y—x) 



X n 



£=±l,n=0 



[d+i(ut-x-y+e{2n+l)L)]^ 



£=±l,n=0 



[5+i{ut-2x+e{2n+l)L)] [S+i(ut-2y+e{2n+l)L)] 



[S+i{ut+x+y+e(2n+l)L)] ^ 

[S+i{ut+2x+£{2n+l)L)] [S+i{ut+2y+£{2n+l)L}] 



and appropriate expressions at finite-temperature. The tunneling DOS is the Fourier trans- 
form of this correlator. 



VII. CONCLUSION. 

We have extended the bosonization technique to a LL connected to resistances computing 
correlators of the boson fields in so doing. The latter are the building blocks allowing 
calculation of the fermion correlators. As side results we derived also the finite-frequency 
conductivity and found that the finite-size LL with IBC is equivalent to an infinite LL by 
virtue of identity of Green's functions whenever impedance matching is realized. 

We also recovered explicitly the properties of an open LL: it corresponds to IBC with 
infinite resistances. 

But in general the LL with IBC has distinctly different properties (it is a dissipative sys- 
tem) and forms a universality class in its own right much as the LL with OBC which exhibits 
critical exponents different from those of the infinite LL. Further interesting developments 

28 



using the results in this paper would be a study of the single particle spectral density which is 
the object of interest in tunneling experiments. A study of the shot noise using the Keldysh 
technique would also be straightforward given the knowledge of the Green's functions. 



Appendix A: Conductivity. 

We use linear response theory: let us consider the perturbation V = — J dxp U where U 
is a voltage. Integrating by parts one gets: V = — J dx E By linear response: 

S {4>{x, t)) = ^ I dy I dt' Gr{x, t- y, t') E{y, t'). 
But ie{t) = ^ therefore: 

x,uj) = --uj I dyGii{x,y,uj)E{y,uj) 







TT TT 



SO that the non-local conductivity is: 

g 

a{x,y,uj) = -ujGR{x,y,uj). 
n 

The expression of Gr{x, y, u) was computed in section lllli This expression allows mea- 
surement of the complex boundary impedances. 

Appendix B: Recovering the Luttinger Liquid with open boundary conditions. 

In this Appendix we will show that the LL with OBC is a special limit of the LL with 
IBC when the reflection parameters at the boundaries are set to unity. Physically this comes 
about because perfect reflection can be equated with having a hard-wall. 

The OBC for the fermion operator: 

'ip{x) = {x = or L) 



I.e. 



MO) = -i^do) 
ML) = -ML) 
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will be derived explicitly from the IBC (with rs — i — td) (we have dropped the phase 
^i2kFL gjj^gg /j.^ — tVtt/L). This will show the perfect equivalence between the LL with OBC 
and the LL with IBC when rg — ^ — td- But already we observe that: 

OBC ^ IBC (rs = 1 = m) 

since cancellation of the fermion operator implies that its current is also zero (and partic- 
ularly its k <^ kp harmonics: see for instance the relation p+{x) = p-{—x) found for the 
OBC, third equation below eq.(lO) of Ref . 4) ). 

But for reflection coefficients = 1 = td the IBC is equivalent to stating that p+ = p_ 
at both boundaries (therefore the current / oc p+ — p_ vanishes). The OBC therefore does 
imply the IBC. 

We now prove the converse and establish: 

IBC (rs^l^rD) ^ OBC. 

The proof actually exists already at 90% in Ref. 5 which in order to quantize the LL 
with OBC actually started from the zero current condition on the boson field: in that work 
the fermionic boundary conditions are used to derive the quantization rules on charges. But 
one can take a reverse standpoint: this allows to derive the fermionic boundary conditions 
by demanding that charge be conserved. (We will also in the course of the proof reconcile 
Ref. 4 and Ref. 5 which find a slightly different Green's function [see eq.(29) of Ref. 5 and 
the discussion after it about the difference with eq.(31) of Ref. 4] .) 

A. Derivation of the OBC. 

1. Mode development of the fields. 

We depart in this Appendix from the definition of the boundaries at x — ±L/2 to set 
them at X — 0, L. For refiection coefficients rs — 1 — Td the IBC imply: 

d^(j)_{L,t) = d^(j)+{L,t). 
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Therefore: 

dtMo,t) = -9t0_(o,t), 

dt<j)-{L,t) = -dtML,t). 

This imphes: 

0(O,t) = Co, <f){L,t)=CL (12) 

where the operators Co and Cl do not depend on time. 

For these values of rs and ro the boson theory is dissipationless and will be described 
by a c = 1 conformal field theory. It will be necessary to make an eigenmode development 
of the fields. To do it we follow closely Ref. 5 who treated the quantization of the OBC 
by using as a starting point eq. ( fT2l ) (i.e. by using the IBC!). Since the field obeys the 
standard wave equation and given the boundary conditions ( eq.(fT2l) ) the mode expansion 
for and G must have the form: 

^{x, t) = <Po + ^Qx + ^/Ky ^^4^ (-^a„e-^^""* + za+e^«'-*) (13) 

n>l 

where qn = nn/ L. 

We impose the standard equal time commutation relations for the fields: 

[0(x); m] = 0, [e(x); eiy)] = 0, [0(x); 0(i/)] = tOix - y) 

where 9{x) is the Heaviside step function. The first two commutators [(p{x); (p{y)] = 
0, [Q{x); Q{y)] = imply that all the commutators vanish except [a„; a^] and [Q; Qq]. 
Expanding the third commutator one gets: 

te{x -y) = — [Q- Go] a; + 2z > a„; al\ 

n>l 

and by making use of the expansion: 

.X „ . sin(g„a;) cos(g„?/) 



i9{x — y) = i— + 2i 



L ^ Tin 

n>l 

which can be proved by using eq. (A-1) of Ref. 4 one finds that the only non-zero commu- 
tators are: 

[ttn, a+] = 1, [Q; Go] = i/v^- 
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We note that 0o is a c-number; this is normal since the operator J which is usually its 
conjugate momentum in the periodic LL does not appear in the theory. We can therefore 
remove it altogether since the LL Hamiltonian is invariant under constant shifts of the phase 
field: (j) — 4> ~ 4>o- 

The zero mode Q is as usual the charge added to the system since the number density is 
p - Po = ^^a;0: 

N-No^ ^(0(L)-0(O)), 
= Q. 



2. Picking the right bosonization formula. 

The usual bosonization formula for the fermion operator used in infinite systems 

ipR/L = ;^^= exp {±4> + 6) 

is not the only one possible; other valid vertex operators for a fermion operator are for 
instance: 



i^R = exp^^/^^{(|) + e) , 

_ g-jax+7 

■0L = — exp iy/n {(j) + Q) 

where a, (5 and 7 are real constants. It is easy to check that these constants do not affect 
the bosonization proof given in Section 5: (3 and 7 correspond to a [/(l)^ x U{1)r invariance 
of the free Dirac lagrangian while a only shifts the chemical potential (substitution in the 
Dirac Hamiltonian leads simply to the additional term a{Npi + Nl)). The effect of these 
constants on the physics is slight: the a term leads to additional oscillations in the Green's 
function, which are in a sense trivial because they only correspond to a shift in chemical 
potential; the other constants have no effect. However when one wishes the operators to 
obey specific boundary conditions they will be necessary. 

Indeed this liberty in the choice of the fermion operator is useful: remember that OBC 
can be obtained from a variety of conditions, for instance Dirichlet BC (vanishing of the 
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fermionic wavefunction) or Neumann BC (vanishing of its derivative). Dirichlet BC = 
at a boundary at x — oi x — L reads: 



MO) = -V'l(o), 

A Neumann BC reads dxip — 0. In the low-energy limit k <^ kp- 

= ^ kp (e^^^^^V'ii - e-^'^'M = 

and therefore at x = 0, L: 

MO) = MO), 

ML) = e-'^'^-'^ML). 

Therefore the same vertex operator can not obey both Dirichlet and Neumann BC at the 
same time (at the same location). We now give the correct prescriptions for both situations 
and also mixed ones (Dirichlet at one boundary and Neumann at the other). 

It will be convenient to define in the following discussion the 'primed' operators which 
differ from the usual ones by a shift of the Fermi vector 



M^,t) = -^^^^ expi^/^^{-(f) + e) exp -i-^{x + ^t). 
These operators correspond actually to factoring out the zero modes from the phase fields: 



t/j'j^{x,t) = —= cxpiy/nOo expi7r^(x — —t) expiy/n {(f)' + 0') 
V 2n6 L K 

V'lI^)^) — ~i= expi^/vrOo exp — i7r^(a; + — t) expiy^ (— 0' + 6') 
V27r5 L K 

(where the prime ' means fields from which the zero modes have been subtracted) . These re- 
lations follow immediately from Campbell-Haussdorf formula and the commutator [Q] 0o] — 
ij^pK {(j)' and 6' commute with the zero modes). 
Dirichlet boundary conditions. 
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In the free-fermion limit the zeros at the two boundaries imply that the Fermi vector has 
the form: kpL = nN so we can remove the phase e~*^^^^ from the BC. 
Let us consider the operators: 



(the upper index D stands for Dirichlet). 

We now use the mode expansion on the fields (eq.(fT3l) ) which implies: 

0'(O) = O = 0'(L). 

This in turn implies: 

^p^{0,t) = —= expi^^Qo exp —in^—t expi^Q' 
Therefore at x = we have trivially a Dirichlet BC: 

^^(o,t) = -v^f(o,t). 

ki X = L: 

ipji{L,t) = -n= expiy/nQQ expmQ expiiT^—{—ut) expiy/nQ', 
ijj^^L.t) = ^= expiy^Qo exp — zttQ expin^—{—ut) expi\^Q'. 

The two expressions are almost identical except for the second term at the right of the 
sign equal. We now enforce charge quantization: the operator Q must have only integral 
eigenvalues. Therefore we have the operator equality: 

QxpmQ = exp —mQ. 

This then implies Dirichlet BC at x = L (with kpL = tcN ): 

^^(L,t) = -^f(L,t). 
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Furthermore using the mode expansion one finds that: 



cj)'{x) + e'{x) 



(l)'{-x) + e'{-x). 



Since: 




1 



expiy/TiQo expiTT— (x — ■ 

Ij 1 



ut 



) exp20F((/.' + e') 




1 



K 



ut 



) expiV^(-0' + e') 



K 



it follows that: 




One can check that the prescription is exactly that of F-G in Ref. 4 (see eq. (9) and eq. 



We have thus recovered Dirichlet boundary conditions starting from the IBC at rs — i — 

vd- 

Neumann boundary conditions. 

Again in the free-fermion limit the zeros at the two boundaries imply that the Fermi 
vector has the form: kpL = nN so we can remove the phase e~*^*^^^ from the BC. 

Since we are at liberty to add constant phases and still get vertex operators for fermions, 
the previous discussion suggests that suitable expressions are: 



Mixed boundary conditions. 

For mixed boundary conditions kp is quantized as kpL = 7r(A^+ |) (which is the correct 
quantization for free fermions in a box with mixed conditions ^ = at one end and dxi^ = 
at the other one). 



(7)). 



This leads immediately to Neumann BC at x — and x — L: 



<(0) 
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- for Dirichlet BC at x = and Neumann BC at x = L: 



MO) = -MO) 

ML) = e-'^'-'^ML) = -ML). 
therefore the expressions of the fermion fields used for Dirichlet BC still work: 

Ki^.t) = ijj'j, 
^PEix,t) = -^^ 

- for Neumann BC at x = and Dirichlet BC at x = L, the Neumann prescription above 
is the correct one. 

A 'Twisted' boundary condition. 

The interested reader may inquire what boundary conditions the standard operators i/jr/l 
can describe; it can be checked that: 

^^/^ = (^'^ + ®) 

V ZTCO 

implies: 

MO) = MO), 
ML) = - ML). 

This can describe a system with any combination of Dirichlet or Neumann BC provided we 
add the condition that one boundary adds a vr phase; a way to do that with free fermions 
is to add a boundary interaction with reflection coeflicient r = — 1 such as 'ip^i'R + i'L^'L — 
ipRipL — ■ If unfold the non-chiral system of length L into a chiral system of length 

2L this corresponds to the theory of a single chiral fermion on a circle threaded by a flux n. 



3. Fabrizio-GogoUn Bosonization of OBC versus Mattsson et coll. Bosonization. 

The earliest theory on the LL with OBC is due to Fabrizio and Gogolin (F-G)^. A 
little bit later the problem was also treated by Mattsson and collaborators (M-E-J)^ with a 
different bosonization scheme namely the standard prescription used for the inflnite system: 

i'R/L = ^^y= exp (±</) + 9) . 
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The treatments lead mostly to the same results although there are in details some minor 
differences. For instance, M-E-J find that the Green's function has an additional modulation 
coming from zero modes not present in F-G's results (see eq. (29) of Ref. 5 and the discussion 
which follows and compare to eq. (31) of Ref. 4): 

Gm-e-j{x, t; y, 0) = e-^^(2no-i)(.-,+.ei^c-^t)G!^_^(^^ f, y, 0) 

(the Green's function written is that of the left fermion ipL with Uq defined by M-E-J as 
no = kpL/T: modi). M-E-,J comment that this difference with F-G will have implications 
for time correlations but do not explain the origin of the discrepancy. 

The previous discussion should hint at the explanation: the additional phase as compared 
with F-G comes from the fact that F-G use the primed operators 

This leads naturally to the question: which is the correct prescription since M-E-J also 
use Dirichlet BC? We now show that M-E-J choice leads to inconsistencies with regard to 
conserved charges and that the correct prescription for Dirichlet BC is indeed that of F-G. 

M-E-J find (eq.(22a-b) of Ref. 5 and eq.(5) and (18) for definitions) for the total charge 

Qc — Qi + Qj.and Qs — Qi — Qj. (after proper rescaling to extract the physical charges): 

2kFL 



Qc = n+l + 



TT 



Qs = m 

where n and m are integers having the same parity. M-E-J derived these constraints by 
imposing the Dirichlet BC. 

But for free fermions with Dirichlet BC at both boundaries: kpL = ttN (which implies 
that ^^-^ is an even number!). M-E-J equations therefore imply that Qc and Qs have 
opposite parity. 

However since Qc — Qs + 2(5 j., Qc and Qs must have the same parity: we have therefore 
a contradiction. Having the correct quantization conditions plays an important role for the 
partition function and for the finite temperature Green's function for the zero modes part 
so the issue is not innocuous. 

An obvious way to cure the problem would be to artificially prescribe that the Fermi 
wavevector is k'pL = t:{N + |) since this avoids the problem with charges: it also removes 
the additional phase M-E-J find in the Green's function and yields Dirichlet BC but of 
course the value of k'p is incorrect in the limit of free fermions. A possible interpretation of 
the shift might be that it proceeds from a change in Maslov-Morse index in the trajectory. 
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B. Recovering the correlators of OBC. 



We will check here directly using the expressions computed with IBC that we recover the 
correlators of OBC. (NB: In order to compare our results with those of the litterature which 
work with x G [0, L], since we worked with x G [— L/2, L/2] we shift all the space arguments 
by L/2 when using the expressions derived in the bulk of this paper.) 

1. Boson correlator. 
Chiral boson correlator: 

For rs — ^ — rD the chiral boson fields correlators read: 



K 



F++(x, t; y, 0) = t)My, 0)) = -4^ E (^^r^)'"' In [5 + i(ut - 2nL - x + y)] . 

n=—oo 

One can also check directly that: 

(0+(x,i)(/)_(y,O)) =- F++{x,t; -y,0), 

{(t>4x,t)(t>+{y,0)) =- F++{-x,t; y,0), 
{<P4x,t)<P_{y,0)) = F++{-x,t; -y,0). 

Using the product expansion of the sine function: 

sin(7r2;) = ttz - — ) 

one gets: 



n>l 



AF++ = t- y, 0) - F++(0, 0; 0, 0) = In 



2LS + liut — X + y) . Triut — x + y) 

— ^ ^ sm ^ — - 

TTO ut — X + y 2L 



This yields: 



i2L ■ n{ut-x+y) 
■kS ''"^ 2L 



for ut — X + y ^ 0; 
for ut — X + y — 0. 



This agrees with the expressions found in the litterature (e.g. eq. (28) in Ref. 5) in the 
limit of zero temperature for finite-length. 
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One checks easily that 

= ^ {^vr^^l^ -S{x-y- ut) + 5(0) 
where S{z) introduced by F-G is (see eq.(A-l) of Ref. 4): 

S{z) = \n—P{z)+if{z) 

TTO 

P z) = - 



2L^/sinh2|| + sin2g 
sinf 



f{z) = arctan ■ , 

exp ^ — cos^ 

This cumbersome expression will be useful to neatly separate the phase in the fermion 
Green's function. 

Non-chiral correlator: 

Using our previous results (eq. (ITTil ) one has (we keep the calculations at zero temper- 
ature for simplicity): 

+00 

(0(x,t)0(2/,O)) = -— J2 In {[6 + z{ut + 2nL)f + {x - yf] 

n=— 00 

K 

+ — \n[6 + i{ut-x-y + {2n + l)L)] 

n=— 00 

K 



+ — ^ \n[6 + i{ut + x + y+{2n + l)L)] 



We regularize by adding — | (0(x, t)^) — | (0(y, 0)^) and by using the infinite product 

2 2 

expressions for cos and sin (cos(|z) = nn>i(l ~ (2n+i)2 ) sin(7rz) = 7rz]^^>^(l — j^) ) 
one gets: 

(0(x, t)0(2/, 0)) - |(0(x, t)2) - 1 (0(y, 0)2) 
= 1^ { (1^)' sin (t (nt -x + y)) sin (^(nt + a; - y)) } 

f / cos(^{nt+x+j/)) cos(^(«t-a;-s/)) \ 1 

47r "^IV cos(fx)cos(fs,) Jj- 

Ref.- does not give explicitly the full correlator for the field and only the following chiral 



correlator for the left field 0l is derived for K = 1: 

pM-E-J ^ ^^^(^^i)^^(^^o))-^(0^(x,t)2>-i(0^(y,O)2) 

= ^HWJ sin(-(nt + a;-,) 

39 



For the sake of comparison with our expression let us rebuild the full correlator using the 
previous equation. Indeed taking into account the OBC leads to: 



F++{x,t;y,0) = t; -y, 0) 

F+_(x,t;y,0) = - F^^{-x,t;y,0) 
F_+{x,t;y,0) = - t; -y, 0) 

which yields finally (after rescaling the fields to make the LL parameteri^' appear): 

F''-^-^ = my, 0)) - I m - \ {<Piy, 0)^) 

= { (1^)' sin (^(nt -x + y)) sin {^{ut + x - y))] 

This is identical with our result if one takes care to shift the origin to the left boundary 
X — s> X + L/2 (since Ref.-^ uses the left boundary as origin). 

Our results are therefore in perfect agreement with the calculations of Ref. 4 and Ref. 5. 



2. Fermion correlator. 

After recovering the boson correlator we turn to the fermion operator for a spinful LL 
(with customary definitions): 

i^R] = ^= expi0F(0| + 61) expi^(x- ^ ^t) 



where the phase comes from the fact that the operators obeying the Dirichlet boundary 
conditions have the zero modes extracted from the exponential as explained above. 
Defining the correlator as: 

Gr] = -i {iljR^{x,t)ij~l^^{y,0)) 
and inserting the expressions of the (charge and spin) chiral fields one finds: 

_ i . TT \ ^ Ui, IT X - 

= "2;^'^P^2L^'^ - ^ - ^^P 2 C'.(x,t; y) 

l/ = C,S u=c,s 
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where: 



C.{x,t; y) = {l-rK-'f^Fl^{x,t; y) {I - K^f ^Fl^{-x,t; -y) 
- y) + Fl^{x,t- -y) 

-\ [F^ix, 0; x) + Fl^{-x, 0; -x) + 0; 2/) + 0; -2/)] } 



Substitution of the chiral boson correlator yields as a function of P{z) and f{z): 

Gm = i^A-A n ^(^ - ^ - ^^¥''P{^ - y + 

^ ' v=c,s 

'P{x + y- u,t)P{x + y + u^ty 



P(2x)P(2y) 
where if we define = exp 2xu- 



expi$(a;, y, t) 



c„ = cosh Xu, St,^ sinh Xi. 
and where the phase factor ^{x,y,t) (it is not the boson field 0!) is given by: 

^x,y,t) = |(x-y)-^^^t + ^J][c^/(x-y-ii,t) 

- slf{x -y + u^t) + s^Cy {f{x + y + u^t) - f{x + y- u^t))] . 

Our expression coincides with that found by F-G (eq. (31) in Ref. 4) apart from two 
things: (i) an unessential constant " T,„k^+k^ which comes from a difference of normal- 
ization of the Fermion operator (eq. (9) of Ref. 4 has a prefactor oc 1/ a/2L while we have 
oc I/VStt^); (ii) the phase factors $(x, y, t) are identical except for the term ^ (F-G 
have ^ J2u it^^ [factor 2 instead of 4] ). 

One can check however by taking the non-interacting limit that the term ^ Ylv k^^ 
correct one: in the non-interacting case the phase must be a chiral function of x — y — vpt 
since we compute a right fermion correlator; substituting Ky = 1, s,^ = and = I , 
Uu — vp in our expression then yields ^{x, y,t) — j{x — y — Vpt) + f{x — y — Vpt) which has 
the right dependence (which F-G can not have with the ^ J^i^ 'k'^ term). 
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In conclusion we have recovered using IBC the fermion correlator for OBC (correcting in 
passing a misprint in Ref. 4). 
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